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I review the results from lattice gauge theory for the properties of the light 1 h exotic 
state and ++ glueball. 

1. Introduction 

QCD is a simple elegant theory, but it is hard to solve it using anything other than 
perturbation theory. Quantities, such as masses, depend on the coupling (g) like M ~ 
e" 1 ' 9 [0], hence perturbation theory can't be used to compute the masses of hadrons such 
as the proton. Our inability to solve QCD non-perturbatively makes it hard to determine 
basic parameters of the quark sector of the standard model. For example, the allowed 
range on the strange quark mass in the particle data table @ is 75 to 170 MeV; a range 
of almost 100%. 

A particularly good test of our understanding of the non-perturbative aspects of QCD 
is to study particles where the gauge field is excited somehow, and hence playing a more 
important dynamic role than in "standard" hadrons. Examples of such particles are 
glueballs (particles made out of the gauge fields) and hybrid mesons (qq and excited 
glue). A very good overview of the interesting issues in hadronic physics, that contrasts 
the hadron spectroscopy approach to the study of confinement with the results from DIS 
type studies, is the white paper by Capstick et al. ||. The only technique that offers 
any prospect of computing masses and matrix elements non-perturbatively, from first 
principles, is lattice QCD. 

I review the results from the lattice for the lightest ++ glueball and 1 h exotic, as 
these states are the closest to being experimentally confirmed. Other recent reviews 0-0 
of lattice results for hybrids and glueballs focus on different aspects of the subject. 

2. Lattice QCD calculations 

Many bound state properties of QCD can be determined from the path integral 

c(t) ~ J dU J dtj) J dip J^O&O^&tye-^-^ (1) 

where Sf is the fermion action (some lattice version of the continuum Dirac action) and 
Sq is the pure gauge action. The path integral in eq. [I] is put on the computer using a 
clever finite difference formalism ]l|] , due to Wilson, that maintains gauge invariance. The 
path integral in eq. p] is evaluated using algorithms that are generalisations of the Monte 
Carlo methods used to compute low dimensional integrals. The physical picture for eq. [I] 
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is that a hadron is created at time 0, from where it propagates to the time t, where it is 
destroyed. The physics from the calculation is extracted using a fit model [[[]]: 

c(t) = aoexp(—m t) + aiexp(-mit) + • • • (2) 

where m (mi) is the ground (first excited) state mass and the dots represent higher 
excitations. Although in principle excited state masses can be extracted from a multiple 
exponential fit, in practice this is a numerically non-trivial task, because of the noise in 
the data from the calculation. Any gauge invariant combination of quark fields and gauge 
links can be used as interpolating operators (0(x,t)) in eq. [|. 

The fermion integration can be done exactly in eq. [L] to produce the fermion deter- 
minant. The determinant describes the dynamics of the sea quarks. In quenched QCD 
calculations, the determinant is set to a constant. Quenched calculations are roughly 
1000 times cheaper computationally than the calculations that include the dynamics of 
the sea quarks. Quenched QCD gives quite a reasonable description of experiment. For 
example, the most accurate quenched calculation of the hadron spectrum, to date, has 
been completed by the CP-PACS collaboration f7|. From the masses of 11 light hadrons, 
they conclude that the quenched approximation disagrees with experiment by at most 
11%. 

In an individual lattice calculation there are errors from the finite size of the lattice 
spacing and the finite lattice volume. State of the art lattice calculations in the quenched 
theory, run at a number of different lattice spacings and physical volumes and extrapolate 
the results to the continuum and infinite volume limit. The increased computational 
costs of unquenched calculations means that most calculations are currently done at fixed 
lattice spacings. One of the most interesting unquenched calculations is being performed 
by the MILC collaboration ||. MILC's calculations include 2+1 flavours of sea quarks 
with a lattice spacing of 0.13 fm, box size of 2.6 fm, and the lightest ratio of the pseu- 
doscalar to vector mass is 0.4. 



3. Results for glueballs in quenched QCD 

Interpolating operators for glueballs are constructed for eq. [I] from closed loops of gauge 
links with the required J PC quantum numbers Some highlights of the results are that the 
lightest glueball is the ++ state with a mass of 1.64(4) GeV ||. The next lightest glueball 
is 2 ++ . The ratio of the tensor to scalar glueball mass is M 2 ++/M ++ = 1.42(6) ||. The 
spectrum of glueball states for other J PC quantum numbers with masses under 4 GeV 
has been comprehensively mapped out by Morningstar and Peardon [p70f] . 

In the real world glueballs will decay to two mesons, hence they will have a decay width. 
Lattice QCD calculations are performed in Euclidean space, for convergence of the path 
integral in eq. |1[ The Euclidean nature of lattice calculations makes the computation of 
inherently complex quantities such as decay widths more involved |TT] . 



The GF11 lattice group computed the decay widths for the decay of the ++ glueball 
to two pseudoscalars |12| to be 108(28) MeV. Although the error bar is only statistical, 
it is encouraging that the width was small relative to the mass, so the ++ glueball may 



exist as a well defined state. The decay widths for individual meson pairs |13| did not 
agree with the predictions from the "flavour democratic" assumption. 
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The experimental situation || for light ++ scalars is very interesting, because there are 
too many states to put into SU(3) nonets, as other particles with different J PC quantum 
numbers, such as the pseudoscalars, can be. The hadrons: / (1370) ,/ (1500) , and 
/ (1710), have masses close to the mass of the quenched scalar ++ glueball. 

In full QCD interpolating operators with ++ can be constructed from quarks and ante- 
quarks, such as ipijj. In full QCD, the pure glue ++ operators will mix with the fermionic 
++ operators. If the mixing is very strong, then the final ++ masses will have little to 
do with the quenched glueball masses. 

Weingarten and Lee |L4| studied the effect of mixing between the glueball and ipip states 



in quenched QCD. They measured the correlation between the ++ glueball states and 
ipip states in eq. [I]. The results were expressed as a mixing matrix 



m g 
E(s) 




(3) 



where m g is the glueball mass, m a (s) is the mass of the non-singlet ++ state at strange, 
and E(s) is the mixing energy. Weingarten and Lee measured: m g = 1648(58) MeV, 
m a (s) = 1322(42) MeV, and E(s) = 61(58) MeV. The qualitative picture that emerges 
is that the / (1710) is "mostly" ++ glueball, and the / (1500) is "mostly" ss. It is not 
clear whether / (1500) being ss is consistent with its decay width |15] . 



The mixing energy E(s) has large lattice spacing errors. For example at a lattice spacing 
of a -1 ~ 1.2 GeV, the Weingarten and Lee WM result is E(s) ~ 0.36 GeV. This has been 



checked by another group's result [[UJ of E(s) ~ 0.44 GeV. 

The analysis of Weingarten and Lee [TJj] depends on the ++ states being well defined 
in quenched QCD. Eichten |TJJ et al. have shown that there is a problem with the non- 
singlet ++ correlator in quenched QCD. The problem can be understood using quenched 
chiral perturbation theory. The non-singlet ++ propagator contains an intermediate 
state of f]' — 7i. The removal of fermion loops in quenched QCD has a big effect on 
the Tj' propagator. The result is that a ghost state contributes to the scalar correlator, 
that makes the expression in eq. |2| inappropriate to extract masses from the calculation. 
Eichten et al. [HJ predict that the ghost state will make the a mass increase as the 
quark mass is reduced below a certain point. This behaviour was observed by Weingarten 
and Lee |I4|] for small box sizes (L < 1.6 fm) for quark masses below strange. It is not 



clear how the problem with the non-singlet ++ correlator in the quenched approximation 



effects the results of Weingarten and Lee |JT4 j , however their most important results come 



from masses above the strange quark mass where the ghost diagram will make a smaller 
contribution that may be negligible 

4. Results for glueball masses in two flavour QCD 

The Weingarten and Lee |14| analysis predicted that the mixing of the ++ glueball 



and ipip states is small. Parts of their calculation have been criticised in JL6|], however, the 



problems with the non-singlet ++ correlator [17| in the quenched QCD will make further 
progress in mixing in the quenched QCD difficult. 

A lattice QCD calculation that included the dynamics of the sea quarks should just 
reproduce the physical spectrum of ++ states. Some insight into the composition of 
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individual ++ states, such as whether a physical particle couples to iftif) or pure glue 
operators, could be studied by looking at the effect of decreasing the sea quark mass. For 
very heavy sea quark masses the theory is more like quenched QCD, where glueballs are 
distinct from i/ji/j operators. 

Hart and Teper |18[ found that the ratio of the ++ glueball mass in rif = 2 QCD to 
the quenched QCD result was: M^ f=2 /M^ uenched = 0.84 ±0.03 at a fixed lattice spacing of 
0.1 fm. The nf = 2 results |18| for the mass of the 2 ++ were consistent with the quenched 



value. As the lattice spacing dependence of the mass of the singlet ++ state in two flavour 
QCD and quenched QCD could be different, a definitive result will only come after a 
continuum extrapolation of the unquenched masses. In quenched QCD 0], the difference 
between the continuum extrapolated mass of ++ glueball mass and the mass at 0.1 fm 
is of the order of 200 MeV. This is the same magnitude of the mass splittings between 
the masses of the experimentally observed particles / (1500) and / (1710). Although the 
current results for singlet ++ states are starting to be interesting, the lattice spacing 
used in unquenched calculations must be reduced before direct contact can be made to 
phenomenology. 

The mass of the ++ glueball on the UKQCD data set are just degenerate with the 



mass of two pions [18]. As the mass of the sea quarks is reduced, two pion states may 



effect the physics of singlet ++ states. 

5. Results for light 1 + exotic mesons 

The quark model predicts the charge conjugation (C = (— 1) L+S ) and parity (P= 
( — 1) L+1 ) of a meson with spin S and orbital angular momentum L. States with quantum 
numbers not predicted by the quark model, such as: J^ x ° otic = are 

known as exotics ||19|| . Exotic states are allowed by QCD. Morningstar and Peardon [10 
claim that there are no glueballs with exotic quantum numbers with masses less than 4 
GeV. 

There are a number of different possibilities for the structure of an exotic state. An 
exotic signal could be: a hybrid meson, that is a quark and anti-quark with excited glue, 
or bound state of two quarks and two anti-quarks (ipipTpip) . 

One possible interpolating operator p0[ , that can be used in eq. [1], for a hybrid 1~ + 
particle is 

0i-+ (x, t) = ip(x, t)"/jF i:j (x, t)ip(x, t) (4) 

where F is the QCD field strength tensor. If F is removed from eq. |4], the operator creates 
the p particle. In this formalism a gauge invariant interpolating operator, for any possible 
exotic hybrid particle or four particle state can be constructed. The dynamics then 
determines whether the resulting state has a narrow decay width, hence it can be detected 
experimentally. In the large N c (number of colours) limit 0,^l|] both exotic hybrid mesons 
and non-exotic mesons have widths that are small compared to their masses. 

There have not been many new calculations of the mass of the light 1 _+ hybrid recently. 
All the results from the various lattice QCD calculations, by UKQCD |22],23|, MILC [g0 



B4fl and SESAM [25] are essentially consistent with the mass of the 1 state around 



1.9(2) GeV H]. The interpolating operators used to create the exotic meson states in 
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the MILC calculations f20j are different to those used in the UKQCD |2| and SESAM 



simulations [B3[] , hence giving confidence that the systematic errors are under control. The 



results for the hybrid masses reported by Lacock and Schilling ||25|| , include some effects 
from dynamical sea quarks. The recent results for the 1 h mass from calculations that 
used an asymmetric @ lattice in time are consistent with the older results. 

There are a number of experimental candidates for light 1 _+ states M. The E852 
collaboration have reported |27j a signal for 1 h state around 1.6 GeV. There is also an 
experimental signal for a 1 _+ state at 1.4 GeV B. 



There has been some recent work |28j on the quark mass dependence of the 1 states. 
The lattice calculations are usually done at large quark masses and the results extrapolated 
to the physical quark masses. The conclusion of J28[ was that the inclusion of the decay 



of the hybrid in the quark mass dependence of the exotic mass could reduce the final 
answer by 100 MeV. The predictions in |28| will be tested as the quark masses used in 
lattice calculations are reduced. The 1~ + state at 1.4 GeV seems low relative to the lattice 
results. 

It is possible that the states seen experimentally are really ipipipip states, in which 
case the operators used in the lattice simulations (eq. f|) might not couple strongly to 
them. Alford and Jaffe |29j studied iptpipip operators with J PC = ++ in a recent lattice 



calculation. The motivation was to gain insight into states such as the /o(980) that some 
people believe is not a ipip meson, but a ipipipip state. A similar lattice calculation could 
in principle be done for the J PC = 1~ + exotic. 

To definitely identify a particle requires both the calculation of the mass as well as the 
decay widths. There has been very little work on strong decays on the lattice. The most 
obvious hadronic process to study using lattice gauge theory is the p — > nn decay, however 



there have only been a few attempts to calculate the g p7T7T coupling [p0 , pl |. Michael 
discusses the problems with the formalism for hadronic decays on the lattice W|. 

In the static quark limit the exotic states on the lattice are described by adiabatic 
potentials. The ground state of the static potential (A lg ) is the familiar Coulomb plus 
linear potential. The excited potential (E u ) is a very flat potential, that can be used with 



Schrodinger's equation to predict the spectrum of heavy-heavy hybrids [§[]. UKQCD [32 
have investigated the de-excitation of the E u potential to the Ai g potential by the emission 
of a light quark loop. In the real world, the decays would correspond to 1~ + — ■> XbV and 
1~ + — > XbS with S a scalar and r\ a pseudoscalar. The decay width of 1~ + — > XbV 
and 1~ + — > XbS transitions were less than 1 MeV and around 80 MeV respectively. The 
various approximations in the static limit mean that these widths have no direct relevance 
to experiment. 

The MILC collaboration |2(| have investigated the mixing between the operator in eq. |] 
and the operator (it <g) ai) eq. [5[ 

^ a i5^y b i5ii^ h (5) 

that has the quantum numbers 1~ + . This type of correlator is part of the calculation 
required to compute the decay width of the 1~ + state to p, and a\. The more complicated 
part is to use eq. ^| in eq. [I] requires some clever numerical work. 
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6. Conclusions 

The glueball spectrum from quenched QCD is essentially complete. The key issues now 
are to quantify the mixing of the glueball states with ipip operators and to determine 
decay widths. 
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